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The motion of self-propelled particles is modeled as a persistent random walk. An analytical 
framework is developed that allows the derivation of exact expressions for the time evolution of 
arbitrary moments of the persistent walk’s displacement. It is shown that the interplay of step 
length and turning angle distributions and self-propulsion produces various signs of anomalous 
diffusion at short time scales and asymptotically a normal diffusion behavior with a broad range 
of diffusion coefficients. The crossover from the anomalous short time behavior to the asymptotic 
diffusion regime is studied and the parameter dependencies of the crossover time are discussed. 
Higher moments of the displacement distribution are calculated and analytical expressions for the 
time evolution of the skewness and the kurtosis of the distribution are presented. 
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I. INTRODUCTION 

Self-propelled particles undergo active Brownian mo¬ 
tion by consuming energy, obtained either from internal 
or external sources. Examples range from the transport 
of motor proteins on cytoskeletal filaments [1] which is a 
biologically relevant system, to the motion of self-motile 
colloidal particles [2] as a nonliving realization. The par¬ 
ticles are powered by the hydrolysis of ATP in the for¬ 
mer case, whereas they use a chemical reaction catalyzed 
on their surface to swim in the latter example. The di¬ 
rected propulsion subject to fluctuations has been de¬ 
scribed by persistent random walks [2-4], where a ten¬ 
dency to move along the previous direction is taken into 
account. A strong self-propulsion overcomes the stochas¬ 
tic fluctuations and directs the motion, which renders it 
ballistic for short time scales [2, 4]. Nevertheless, the 
interplay between self-propulsion and random motion in 
general may lead to various scenarios of anomalous dif¬ 
fusive dynamics on varying time scales [4]. The influence 
of self-propulsion diminishes over time and eventually a 
crossover to an asymptotic diffusive regime occurs. 

Even in the absence of self-propulsion, the stochastic 
motion of particle may remain complicated because in 
general a random walker can perform steps with arbi¬ 
trary turning angles and variable step lengths. More¬ 
over, there can be a relation between the step size and 
the turning angle of each step. Generalized random walks 
had been studied e.g. in the context of animal and cell 
movements as a Markovian process [5, 6], i.e. by consid¬ 
ering the motion as a series of independent draws from 
the step-length and turning-angle distributions for each 
step. While the focus of prior studies has been more 
on the asymptotic diffusion coefficient of such random 
walks, the short-time behavior is neither thoroughly in¬ 
vestigated nor completely understood. 
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In the persistent random walk model that we study 
here a particle moves straight in continuous space in a 
randomly chosen direction over a randomly chosen dis¬ 
tance and then changes direction by a randomly chosen 
turning angle. The typical trajectories of the random 
walker depend strongly on the chosen turning angle and 
distance distributions. For small values of the angular 
change the new direction will be strongly correlated with 
the old direction, introducing a directional memory into 
the model without changing the Markov property of the 
process. The emerging intermittent directional bias is 
controlled by the characteristics (mean, width, asymme¬ 
try, etc.) of the turning angle distribution and the prob¬ 
ability with which the direction is unchanged, i.e. the 
processivity. The bias decays with time after a few turns 
and the directions of the particle motion become asymp¬ 
totically randomized. 

Memory effects have also been included in other ran¬ 
dom walk models. In fractional Brownian motion [7], 
sub or superdiffusive motion is observed asymptotically. 
In continuous time random walk (CTRW) models [8], 
for example, true non-Markovian effects can be imple¬ 
mented via broad waiting time distributions which lead 
to a subdiffusive dynamics of the walker. Although these 
approaches are conceptually very exciting, a direct com¬ 
parison to experimental results is sometimes difficult. In 
case of intracellular transport, for example, spatial con¬ 
finement as well as a finite observation time imply that 
the asymptotic behavior of the walker is not accessible. 

Here, we develop a general analytical framework to 
study persistent random walks over the whole range of 
time scales. The goal is to clarify and disentangle the 
combined effects of self-propulsion p and the stepping 
strategy of the walker, consisting of its step-length 
and turning-angle R((f>) distributions. The method en¬ 
ables us to analytically determine the time evolution of 
arbitrary moments of displacement. Using this approach, 
the second moment, i.e. the mean square displacement, 
has been recently analyzed [4] revealing a variety of sig¬ 
natures of anomalous diffusion on short time scales even 
in the absence of viscoelasticity, traps, or overcrowding. 



2 


These elements were frequently identified in the nature 
of the biological environments and received considerable 
attention as the possible sources of subdiffusion [9]. 

An alternative way to interpret the model and its out¬ 
come is to consider the motion on complex networks, such 
as motor proteins on the cytoskeleton. These motors 
have an effective processivity p (i.e. a tendency to move 
along the same filament [ 10 - 12 ]) and may switch to a 
new filament at the intersections of the network. From a 
coarse-grained perspective, one can describe such motion 
by a persistent random walk on the nodes of the net¬ 
work, thus, P(£) and R(4>) represent the distributions of 
the segment-length £ between neighboring intersections 
and the angle (j) between intersecting filaments, respec¬ 
tively. Note that the cytoskeleton is a dynamic network 
due to the underlying growth and shrinkage of filaments, 
thus, the structure on which the transport takes place 
often changes. This justifies the relevance of the pro¬ 
posed stochastic approach, where the network structure 
is always implicitly given. 

Our analytical approach provides a recipe to ob¬ 
tain any arbitrary moment of displacement, yet we ex¬ 
tend the calculations to the third and fourth moments, 
which are of particular interest, e.g. in the evaluation 
of the skewness and kurtosis, or to obtain the variance 
tr r 2 =(r 4 ) —(r 2 ) 2 (and thus the standard error) of the 
mean square displacement (MSD) of a persistent ran¬ 
dom walker, ay 2 is also a useful quantity for estimating 
the first moment (r) of the net displacement: In the ab¬ 
sence of an exact analytical solution for (r), an approx¬ 
imate expression (r) ~ \J (r 2 )(l—| ^ 2 ) was proposed 
[13] by means of a Taylor expansion of the square root 
function. The asymptotic behavior is however shown to 
follow (r) = i v ,/ 7 r (r 2 ) [14]. 

We first introduce the general approach to obtain ar¬ 
bitrary moments of displacement. Next, an analytical 
expression is derived for the special case of the mean 
square displacement in two dimensions. This part is the 
full exposition and expansion of the results presented in 
[4]. Then we clarify the similarities and differences be¬ 
tween the results in two and three dimensions for per¬ 
sistent walks which are symmetric around the previous 
direction of motion, and briefly discuss the case of asym¬ 
metric turning-angle distributions in two dimensions. Fi¬ 
nally, the calculations are extended to higher moments 
and cumulants of displacement, and the probability den¬ 
sity of the position of the random walker and its time 
evolution are investigated. The paper is organized in the 
following manner: Section II contains the description of 
the master equation formalism for the persistent motion 
of random walkers in two dimensions. In Sec. Ill, we dis¬ 
cuss symmetric persistent motions in three dimensions. 
The analytical predictions for MSD are compared with 
simulation results in Sec. IV. The calculations are ex¬ 
tended to higher moments and cumulants of displacement 
in Sec. V, and the results are compared with simulations. 
We investigate the parameter dependence and time evo¬ 
lution of the probability distribution of the net distance 


from the origin in Sec. VI. Moreover, the coupling be¬ 
tween longitudinal and perpendicular transport is briefly 
discussed, and an analytical expression is obtained for 
the probability that the direction of motion after n steps 
makes an angle a with the current direction of motion. 
Section VII concludes the paper. 


II. MODEL 


We consider the motion of a self-propelled particle as 
a persistent random walk consisting of steps of different 
lengths and orientations. Here, a two-dimensional mo¬ 
tion is introduced and the extension to three dimensions 
is discussed in the next section. The stochastic motion 
is described in continuous space and discrete time as fol¬ 
lows: At each time step, the particle takes a step of length 
l along either its previous direction with probability p or a 
newly chosen direction with probability s=l—p, as shown 
in Fig. 1. Thus, p represents the self-propulsion of the 
particle. We assume probability distributions R{4>) and 
P(£) for the rotation angle <f>=9— 7 and the step length 
l of the walker, respectively. The following master equa¬ 
tion expresses the evolution of the probability density 
P n (x, y\9) for the particle to arrive at position ( x 1 y ) along 
the direction 9 at time step n: 


P n+ i(x,y\9) =pJd£P(£) P n (x-ecos(9),y-£sm(9)\9) 

+ sJd£P(£) J d"/R(9—'y)P n (x—£cos(9),y—£sin(9)\"/). 

(!) 

The terms on the right hand side of the above equation 
correspond to persistent motion with probability p and 
turning with probability s. One can obtain arbitrary 
moments of displacement since they are accessible by the 
derivatives of the Fourier transform of P n (x , y\9), which 
is defined as 


P n (u)\m) 




P n (x,y\9). 


( 2 ) 


The arbitrary moment ( x kl y k2 ) is given by 


(x kl y k2 ) n = Jd9 jdy Jdx x kl y k2 P n (x,y\9) 


= ( ,) fcl +fc 2 d fcl+fc2 P n {u x ,uj y \m=0) 


duj kl duj k2 


( 3 ) 


(cOx )* J dy') — (0j0) 


To study the diffusive behavior of particles one deals with 
the first and second moments of P n (x,y\6), namely (x), 
(y), (x 2 ) and (y 2 ). Thus, we first focus on the derivation 
of these quantities in the following. The same procedure 
is followed in Sec. V to obtain higher moments of dis¬ 
placement. A similar Fourier-Z-transform technique was 
applied to study diffusive transport of light in foams [15]. 
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FIG. 1. (color online). Trajectory of the walker during two 
successive steps. 


The first two moments along the x direction are given by 


.dP n (ui x ,ui y \m=0) 

(4 = -i --- 


8uj r 


(uJx i^y ) — (OjO) 


(X)n = (-0 


2 d 2 Pn{u x ,Uy\m=0) 


dull 


(4) 


(U! x ,Uy)=( 0,0) 

Similar expressions can be written for ( y) and ( y 2 ). 
Fourier transforming Eq. (1), we find 

OO 

P n+1 (ui,a\m) = E 


k=—c 


i k e~ ika x 


P n (co, a\m+k) (p + s lZ(m+k)) Jd£ £F(£)Jk(oj£) 


(5) 


where J fe (z )=^ J^d0 e izcose e- ikd is the fcth order 
Bessel’s function and 7 Z(m)= fl n d(j) e lrn ^R{(j)) is the 
Fourier transform of the rotation angle distribution R((f>). 
One can expand P n (ui,a\m) as a Taylor series 

P„(w, a\m) = Qo,n(a\m ) + iu> (£) Qi,n(a|m) 

- \u 2 (£ 2 ) Q2,n{ot\m) H-, 

with (£) and (£ 2 ) being the first and second moments of 
the step-length distribution P(£). From Eqs. (4) and (6) 
we obtain the moments in terms of the Taylor expansion 
coefficients Qi^ n (a\m): 


(x) n = fd£F{£)£Q Un ( 0|0) =<*>Q lin (0|0), 

<»>„ = fd£P(£)£Q hn {^-\0) =(£)Q hn (l | 0 ), 

J * (7) 

(x\ = Jd£T(l)t 2 Q 2 ,n(0|0) = (£ 2 )Q 2 ,n (0 10), 

(y\ = Jd£P(£)£ 2 Q 2 ,n{\ |0) = <£ 2 >Q 2 ,n(f |0). 

By Taylor expansion of both sides of Eq.(5), one can 
collect all terms with the same power in ui, leading to the 
following recursion relations for the Taylor coefficients 
Q 0 .n(a\m), Qi >n {a\m), and Q 2 ,„(a|m): 

Qo,n+i{a\m) = Qo,n{a\m)(p + sTZ(m)), (8) 


Qi,n+i(ot\m) = Qi, n {a\m)(p + sTZ{m)) 
1 


e ia Qo,n(ct\m-l)(p + s77(m—1)) 


e ia Qo,n(,a\m+l)(p + slZ{m+l)) 


(9) 


Q 2 ,n+i(a|m)= 

{n: 


■^Qo,n{cx\m) + Q 2 ,n(a\m) (p+sU(m)) 
e ia Qi,n{a\m-l)(p + sTZ(m-l)) 


+ e ia Qi, n (a\m+l)(p + sTZ(m+l)) 
+ \e 2ia Qo,n(a\m-2)(p + sTZ(m-2)) 


( 10 ) 


e 2lot Qa,n{pi\m+2)(p + sTZ(m+2)). 


The coupled linear equations (8), (9) and (10) can be 
solved by means of the ^-transform technique. The z- 
transform G(z) of a function G n of a discrete variable 
n = 0,1,2, ■ ■ ■ is defined as 


OO 

G{z) = Y J G n z~ n . (11) 

n —0 

By applying the z-transform to Eqs.(8)-(10), one obtains 
a set of algebraic equations for Qo(z,a\m), Qi(z,a\m) 
and Q 2 {z,a\m) quantities (see Appendix A). 


A. The quantities of interest 

The main goal is to evaluate (x) n , {y) n , ( x 2 ) n , and 
{y 2 ) n , which are given in terms of Qi, n (o:\m) quantities 
in Eq. (7). Here, we derive (x) n and ( x 2 ) n , and a similar 
approach can be followed to obtain (y) n and ( y 2 ) n . The 
z-transform of Eq. (7) leads to the following expressions 
in 2 -space 

OO 

<*>(*)= E W Ql,n(0|0) = W Ql(z, o|o), 

»° ( 12 ) 

{x 2 )(z)= E z~ n {£ 2 ) Q2,n(0|0)=(^ 2 ) Q 2 (z, 0|0). 

n= 0 

From Eqs.(12), (A2), and (A3), one obtains the first and 
second moments of x in the z-space (See Eqs.(Bl) and 
(B2) in Appendix B). The last step to get the moments 
(. x) and (x 2 ) in real time is the inverse z-transforming of 
the z-space moments [i.e. Eqs.(Bl) and (B2)]. By intro¬ 
ducing Ai=p+sTZ(i), the resulting moments are: 

(x) n = <*>Qi,„=o(0|0) (13) 

m l—A n , 

H—2 _ Qo,n=o(0| — 1 )A_i 14 i 

+^Qo,„=o(0|l)A 1 ^, 
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{x\ 


2_ A n y _ A n 

+ (^) 2 Ql,n=o(0| — 1 )-A _1 ^ ^ + W 2 Ql,n=o(0|l)Ai Tj—^ 

m 2 4 n J 2 + A",(A i-i) + A 2 - a 

+^-Qo,n=o(0| - 2) a_,a_ 2 X -At — r^TJ — l Wl ~ ^ 

2 (A_ 1 -A_ 2 ) (A_ 2-1) (l-i4_i) 

+ ^!n rol2M A :: <4 + ^(A-l)+4-A 
2 g °’"-“ (0|2 >' 4 '' 4 - X (A,-A,) (21,-1) (1-20 

(£) 2 . , .[ A n _!~ A_in + n-l A" — A n + ra— 

+ — Qo,»=o(0|0) (1 _ A _ i)2 + A 1 (1 _ Ai) 2 


A^-A.n + n-l' 


U 2 ) 1— A n o /£ 2 \ 1-A 

A - — Qo,n=o(0| - 2)A~2^_a + Qo,7l=o(0|2)A 2 


(l-AJ 2 

1-A” 

2 


(14) 


B. Isotropic initial condition 


For the isotropic initial condition 
Pq(x, y\0)=T;${x)6(y), one finds from Eq. (2) that 
Po(w|m)=^/^d0e ime =2i^. Then, using the 
expansion equation (6) for w=0, it can be seen that the 
only nonzero Q quantity is Qo,n=o(c> ! |0)=l (for m= 0). 
Therefore, Eq. (14) leads to 

<*>„ = 0, (15) 


and Eq. (14), after replacing s with 1 —p, reads 

(*\ = 

{(■? 


A + 


{y+K(l)- P K{l)) ( P +1Z(-1)- P 1Z(-1)) 

(l-p)(l-K(l)) (i- p )(i-K(-l)) 


(l) 2 * (p+K(l)- P K(l)) 

2 (l-p) 2 (l-P(l)) 2 L 

(l) 2 (p+1Z{-l)-p1Z{-l)) 

2 (l-p) 2 (l-P(-l)) 2 


(p+H(l)-pTZ(l)) n -l 

(p+K(-l)- P K(-l)) n -l 


(16) 


with \=(l 2 ) / (l) 2 being the relative variance of the step- 
length distribution. The y-component of the mean 
square displacement ( y 2 ) n has the same form as shown 
in Eq. (16) due to symmetry. Therefore, one obtains 
{r 2 ) n =(x 2 ) n +(y 2 ) n =2(x 2 ) n . 


C. Long-time behavior 


From Eq. (16) in the limit of long time (i.e. (Hoo) one 
obtains the asymptotic mean square displacement as 


(r 2 )n/(t) 2 * 

\. (p+n i)-p^(i)) 

+ (i-p)(i-n(i)) 
(p+K(l)-pK(l)) 


(p+n{-i)- P n(-i)y 

(p+n(-i)-pTi(-i)) 


(l-p) 2 (l-Tl(l )) 2 (i—p) 2 (i—iz(—i )) 2 ' 


Assuming that the particle moves with a constant speed 
v during the ballistic parts of motion, the elapsed time 
after n steps is r = n(£)/v. The diffusion constant D is 
related to ( r 2 ) n as 

(r 2 ) n = 4 Dt, (18) 

thus, we find 

’ { P +lZ(l)-plZ(l)) (p+K(-l)-pK(-l)) - 
(l-p){l-K(l)) (l-p)(l-ft(-l)) 

(19) 



D. Turning with left-right symmetry 


The distribution R{<j>) reflects to what extent the di¬ 
rections of the successive steps are correlated. The ana¬ 
lytical method presented in this section allows us to han¬ 
dle an arbitrary function R(4>), however, we are particu¬ 
larly interested in the distributions with equal probabili¬ 
ties to turn clockwise or anticlockwise. This implies that 
1Z(1)=1Z(—1)(=1Z). The asymmetry of the turning angle 
with respect to the arrival direction is quantitatively re¬ 
flected in the value of 7Z which ranges between —1 and 
1, with zero denoting a uniform case and negative (posi¬ 
tive) values corresponding to a higher chance of motion to 
the near backward (forward) directions in the next step. 
When left-right symmetry holds the imaginary part of 
lZ(m) vanishes, thus, TZ becomes 

7 Z= f d(j) cos(</>) R(4>) : (20) 


and Eq. (16) reduces to 


{x 2 ) n =\n{l) 2 


A + 


2 (p+lZ—pTZ) 


-w 


(l-p)(l-rc)J 

2 (jp+1Z—pTZ) 


(1 ~P) (1 -K)‘ 


(p+IZ-pTZ) 



(17) 
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III. EXTENSION TO THREE DIMENSIONS 


The analytical approach of Sec. II can be straight¬ 
forwardly generalized to the persistent motion in 
three dimensions by introducing the probability den¬ 
sity P n (x,y, z\(f> : ip) for the particle to arrive at position 
(x, y , z) at time step n along the direction characterized 
by the azimuthal and polar angles <f> and tp, even though 
the calculations for the general motion in three dimen¬ 
sions are quite lengthy. However, the processes with sym¬ 
metric turning-angle distribution are of particular impor¬ 
tance since usually the rotational symmetry holds in bi¬ 
ological applications. Thus, we restrict ourselves in this 
section to the turning-angle distributions with cylindrical 
symmetry with respect to the incoming direction. This 
is the three dimensional analogue to those processes in 
two dimensions which obey left-right symmetry. 

Similar to the 2D case, we introduce R((/>) as the prob¬ 
ability of turning with angle (f> with respect to the in¬ 
coming direction (see Fig. 2). The polar angle tp is sup¬ 
posed to be uniformly distributed over the range [0, 2tt\. 
Note that, in contrast to the 2D case, the normalization 
condition of R((f>) in 3D requires an integration over all 
possibilities of ip i.e. R{4>) sin(</>) dcp=l. Here again the 
corresponding Fourier transform of R((/>) is real. Since ev¬ 
ery two successive steps lie in a plane, as shown in Fig. 2, 
one can intuitively write the same two-dimensional mas¬ 
ter equation [Eq. (1)] to describe the motion in three 
dimensions. By solving this master equation, one gets a 
similar expression for ( x 2 ) as the 2D solution presented 
in Eq. (16) in the case of TZ(1)=TZ{— 1), even though with 
different prefactors: 


{x 2 ) n =n 


(ey 


A + 


2(p+£—p£) 

(M( i-*)J 


2<£} 2 (p+£-p£) 


(p+£-p £) n -1 


( 22 ) 


3 (1 -p) ( 1-0 L 

Here, £ is the real part of the Fourier transform of the 




FIG. 2. (color online). Illustration of symmetric rotations 
with respect to the incoming direction, in two (left) and three 
(right) dimensions. 



FIG. 3. (color online). Time evolution of the mean square 
displacement ( r 2 ) scaled by n, for different values of A, p, and 
TZ. The solid lines correspond to analytical predictions via 
Eq. (21), and the symbols denote simulation results. 


rotation-angle distribution 


£ = f dcj) cos(^) R{(j>) sin((/>). (23) 

Jo 

Finally, one can obtain the total mean square displace¬ 
ment ( r 2 ) n = (x 2 ) n +(y 2 ) n +(z 2 ) n , which has the same 
form as in 2D, only TZ is replaced with £: 


(r 2 ) n =n(£) 2 


-{t ) 2 


\ 2(p+£-p£) 

(i-r)(i-bJ 

2 (jp+£—p£) 


(1-p) (1-£Y 


n (24) 


(p+£-p £) n -1 


IV. SIMULATION RESULTS FOR MSD 

In this section we compare the analytical predictions 
with the results of extensive Monte Carlo simulations ob¬ 
tained from the same step-length T(£) and turning-angle 
R(4>) distributions, and self-propulsion p. The formalism 
introduced in Secs. II and III enables us to handle any 
arbitrary function for R{<j>) and F{£), nevertheless, we 
restrict R(4>) in this section to symmetric distributions 
along the incoming direction for simplicity and because 
of its practical applications in biological systems. 

We first investigate the overall behavior of the mean 
square displacement for different values of A, p, and TZ. 
A is a measure of the heterogeneity of the network struc¬ 
ture or the diversity of the step sizes, and TZ quanti¬ 
fies the anisotropy of the structure or the asymmetry of 
the turning angles of the walker. The characteristics of 
F(t) and R{<f>) distributions can be considered as the 
stepping strategy of the random walker which may be 
tunable externally (e.g. by controlling the external agi¬ 
tation imposed on a driven granular system [16]) or in¬ 
ternally (by controlling the strength, density, and spatial 
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FIG. 4. (color online), (a) The crossover time n c to the 
asymptotic diffusive regime versus self-propulsion p. (b) The 
asymptotic diffusion coefficient D in terms of self-propulsion 
p (light gray curves) or anisotropy 7 Z (dark gray curves). The 
solid, dashed, and dotted lines correspond to A=l,2,5, re¬ 
spectively. 


arrangement of obstacles in the system, or by adjusting 
the underlying structure of the environment such as a 
porous medium [17]). p is the self-propulsion of the par¬ 
ticle (equivalently, the processivity or persistency of the 
walker). The case p=7Z=0 and A=1 corresponds to a 
simple diffusion (see Fig. 3). When p and 1Z are both 
positive, they cooperate to send the walker to the near 
forward directions more frequently, resulting in superdif¬ 
fusion at short time scales. If 7Z is negative, it competes 
against p which may lead to sub, normal, or superdiffu¬ 
sion. At the extreme negative value of 7 Z (i.e. TZ—> — 1), 
an oscillatory phase can be observed where the particle 
experiences a nearly back and forth motion [4]. 

It can be seen from Fig. 3 that the asymptotic behavior 
of all curves is diffusive. This is due to the fact that there 



FIG. 5. (color online). Time evolution of (r 2 ) for two different 
distributions J-(€) with the same moments (£) = 2 and (£ 2 )=5 
(A=1.25). Inset: Comparison between the exponential distri¬ 
bution J r (£)=e 1 ~ e (£ E [l,oo]) and the uniform distribution 

(^ min=0 .268, Anax=3.732). 

v ' *-max l min v ' 



FIG. 6. (color online), (a) Mean square displacement (r 2 ) 
vs n for three different distributions with 1Z=0 given in 
the text. The step-length distribution is chosen to be either 
J-(l)=5(£— 1) (with p=0) or J r (£)=e 1 ~ e (with p=0.9). Insets: 
The possible directions of motion at the next step for each 
R(<t>). (b) A few sample distributions R(4>) with 7?.~0.4 (top) 
and 7?.~—0.4 (bottom). The dotted lines show the arrival 
direction and the arrows indicate the possible directions in the 
next step, with length being proportional to the probability. 


is no preferred direction in the system and the effective 
correlations which exist between successive step angles 
are short-range. The crossover time n c to asymptotic 
diffusion can be estimated by balancing the linear and 
exponential terms in Eq.(21). In Fig. 4(a), n c is shown 
as a function of self-propulsion for several values of 7 Z 
and A. Increasing p and/or 7 Z delays the crossover, while 
the walker gets randomized more quickly for strong het¬ 
erogeneities. The asymptotic diffusion coefficient varies 
by several orders of magnitude with control parameters 
p, 7 Z, and A [see Fig. 4(b)]. 

A remarkable outcome of the analytical formalism is 
that the anomalous diffusive motion of the particle is 
fully described by the self-propulsion, and the character¬ 
istics of the step-length and turning-angle distributions, 
namely the first two moments (£) and (£ 2 ) of J-(£) and 
the Fourier transform of i?(</>). Therefore, one expects 
that stepping with different distributions but with the 
same key characteristics mentioned above should lead 
to the same results, independent of the functional form 
of the distributions. To verify this finding by Simula- 
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ric distribution over the range [— 7 r/ 6 , 7 r/ 3 ] which cre¬ 
ates clockwise spirals. Here 7£(+l)~0.87+i 0.23 and 
TZ(— 1)~0.87—i0.23, thus, 1). The asymp¬ 

totic diffusion coefficient [Eq. (19)] is however a real 
number D/v(£)=\ [A+2 ] in the absence of self¬ 

propulsion, with A and B being the real and imaginary 
parts of 7£(±1). We obtain D/(v(£)/4) ~ 1.0, 43.4, and 
2.7 for R\{(f>), R, 2 {<t>), and i? 3 (</>), respectively, in agree¬ 
ment with the simulation results. 


FIG. 7. (color online). Typical trajectories at p=0 and A=1 
for isotropic (left), forward symmetric (middle), and forward 
asymmetric (right) turning-angle distributions, after the same 
number of steps. The arrows show possible directions of mo¬ 
tion in the next step. 


tions, we first choose an isotropic distribution R((f>)= 
and compare two different step-length distributions with 
the same (£) and (£ 2 ) moments. As shown in Fig. 5, 
the simulation results match remarkably for an expo¬ 
nential function J-{Jt)=e 1 ~ L and a uniform distribution 
J r (£)= i/ ^~ £m i n ' )+g ^ max ~ f) ~ 1 (H(x) is the Heaviside step 

v tmax f-min 

function), both with (£)=2 and (£ 2 }=5. 

Next, we choose a given step-length distribution (ei¬ 
ther P(£)=5(£— 1 ) or P(£)=e 1 ~ e ) and compare three 
different turning-angle distributions: a uniform func¬ 
tion Ri (</>)= 27 , a motion restricted to left or right 
directions R 2 {(f>)=^(S((j)—n/ 2 )+S((f>—Tr/ 2 )), and a mo¬ 
tion restricted to forward or backward directions 
R 3 ((/))=^(5(</))+5(fi — 7r )). All these examples correspond 
to 1Z=0 i.e., on average, they have no preference for for¬ 
ward or backward motion. Figure 6 (a) reveals that there 
is a perfect agreement between the simulation results 
obtained for these different turning-angle distributions. 
One can also generate positive or negative values of 1Z 
from different R{(f >) distributions. Several examples are 
shown in Fig. 6 (b) for 7?.~0.4 (or 7£~—0.4), which all 
lead to the same diffusive motion. 

So far, only symmetric distributions are studied. Now, 
we briefly investigate turning-angle distributions which 
are asymmetric with respect to the incoming direction. 
Let us consider two-dimensional walks for simplicity. 
An asymmetric R((f>) in this case means that the left- 
right symmetry of turning is broken, leading to (anti-) 
clockwise spiral trajectories. A comparison is made in 
Fig. 7 between the trajectories obtained from three dif¬ 
ferent uniform distributions over the range (</> m in, </>max): 
Ri{<f>) is an isotropic function corresponding to a nor¬ 
mal diffusion (1Z=1Z(+1)=1Z(— 1)=0), i? 2 (</>) is a sym¬ 
metric function (</> m in=— 7 r/ 6 , </>max= 7 r/ 6 ) which results 
in TZ=TZ(+1)=TZ{— 1)~0.95, and Rz{(j)) is an asymmet- 


V. HIGHER MOMENTS AND CUMULANTS 


The procedure described in Sec. II enables one to ob¬ 
tain any arbitrary moment of displacement. To better 
clarify the proposed recipe, we extend the calculations to 
the third and fourth moments in this section, which are 
sufficient to derive up to the fourth cumulants of displace¬ 
ment and obtain the skewness and kurtosis of a persistent 
random walk which are measures for the asymmetry and 
peakedness of the probability distribution, respectively. 
We also compare the analytical predictions with Monte 
Carlo simulation results. 

From Eq. (3), the third and fourth moments of the 
displacement are given by 


= ( ;)3 d 3 P n (u J X ,U}y\r n 


=0) 


dul 


(w x ,i i i a ) = (0,0) 


(25) 


and 


= H) 


4 d 4 P n (uj x ,u) y \m=0) 


dut 


(j-^x i^y') — (0 » 0 ) 


(26) 


Moreover, by expanding P n (u), a\m) up to the forth order 
terms in w one finds 


P n (w, a\m) = Qo,„(a|m) + iu> {£} Qi, n (a\m) 

- \u 2 {£ 2 ) <52,n(a|m) 

- 7W 3 {£ 3 ) Q3,n(a|m) 

0 

+ ±. 4 (£ 4 ) Qi, n (a\m) H-, 


(27) 


which results in the following relations between (a; 3 ) or 
(a; 4 ) and the Taylor expansion coefficients 


{x\ = Jd£P (£)£ 3 Q 3i „( 0|0) = (^ 3 )Q 3 ,n(0|0), 
(x\ = Jd£P (£)£ 4 Q 4 ,n(0|0) = (£ 4 )Q 4 ,n(0|0). 


(28) 


Then, by following the procedure introduced in Sec. II 
one obtains recursion relations for < 33 ,ri(ct|m) and 
Q 4 , n (a\m) as 








Q 3 , n +i(a\m) = 


Qi,n(oL\m)+Q3 t n(a\m) 


3 um 

2 (£3) 

(DU 2 ) l 

^ Q 2 ,n(a\m+l) + -Qo,n(a\m+l ) 

(DU 2 ) 1 

^ Q2,n(oi|m-l) + -Qo, n (a\m-l) 


( p+sH(m )) 


x 2 e M (p+s72-(?n+l)) 
x ^e l “(p+s7?.(m— 1)) 


+ 


^^0 1 ,„(a|m+2) 


U 3 ) 


Qi,n(ot\m— 2 ) 


e 2ia (jp+sll(m+2)) 
e 2l “ (p+s TZ(m—2)) 


3l “Qo,n(a|m+3)(p+s7?.(m+3)) 


+ -e 3la Q 0)n (o;|m—3) (p+s lZ(m—3 )), 


(29) 


and 


Q 4 ,n+i(a|m) = 


+ e 


3 (£ 2 ) 2 

-Qo,n(a\m)+3^yQ2,n(a\m)+Q4:, n {a\m) 

2^4^Q 3 ,n(a|m+l)+|^^Q lirl (a|m+l) 


+ e 


+ e 


v-’ 4 ) 

3 (t 2 ) 2 


2 (D) 

Q3,n(a\m-l)+^ ( DU ■ ^ Q 1[n (a|m-1) 


2 (£ 4 ) 


2 ^ 4 * ,( 32,ra(Q!|m+2) + -Qo,n(a|77l+2) 

3 (£ 2 ) 2 1 

2 ^- ( ? 2 ,n(a|m- 2 ) + -Q 0in (a|m-2) 


(p+s7?.(m)) 

(p+s 7?.(m+l)) 
(p+s7?.(to- 1)) 
(p+s7£(m+2)) 


(p+s7?.(m—2)) 


1 -3iaM!> 

2 (D) 


i 


<5i in (a|m+3)(p+s7?.(TO+3)) + -e' 


. (DU 3 ) 


3 ia __ 

2 (£ 4 ) 


Qi,n(a|"i—3) (p+s TZ(m—3)) 


+ 4M Qo,n(o:|m+4)(p+s7?.(?n+4)) + ^e ila Q 0 , n (a\m-4:)(p+sll(m-4)). 


(30) 


The corresponding algebraic equations for Q 3 (z, a\m) 
and Q 4 (z,a|m) after the z-transform are given in Ap¬ 
pendix A. From Eq. (28), the third and fourth moments 
of x in the z-space can be obtained as 

OO 

(x 3 )( Z )= J2 Z~ n U 3 ) Q3,n(m = U 3 ) Q*(z, 0|0), 

»° (31) 

(x A )(z)= Y z~ n U 4 ) Q*Am=U 4 ) Qi(z, o|o). 

n —0 

By inserting Q 3 (z, 0|0) and Q 4 (z, 0|0) from Eqs. (A4) and 


(A5), and inverse z-transforming of the z-space moments, 
we finally obtain (x 3 ) n and (x 4 ) n , which are very lengthy 
equations. However, they fortunately reduce to simpler 
forms when we consider the most interesting cases. For 
example, the isotropic initial condition leads to 


<z 3 }„ = 0, (32) 

and if we further limit the motion to a constant step size, 
T(i)=5(l—L ), and turning with left-right symmetry, the 
resulting (x A ) n reads 





































9 


( x \ = J^L 4 


2(A 1 +A 2 ) 2 A 2 +1 4A” + 1 (n+l)^A 1 (3A 1 +l) — ^ 2 (^+ 3 )^ 


(i-A 2 y(A-A 2 y (i - Ai y{a-a 2 ) 

4A™ +1 ( - 8(2+4, )A\A 2 +{A 1 (7A 1 +4)+l) A 2 i +(A 1 (34, +8)+l) 4 2 ) 

( n +l)(l+^i)( - 34 2 + (4,—3)(54 1 +3)4 a +124 1 +7^ 2 (n + l)(n + 2)(l+4,) 2 


(1-A 1 )3(1-A 2 ) (1-4J 2 

2 A,A 2 (A 3 +18A 1 +24)-A 2 (A 1 +l)^A 1 ((A 1 -5)A 1 +23)+5^+A 1 (A 1 +2)((A 1 -6)A 1 -10)+10A 2 -3' 

(l-4J 4 (l-4 2 ) 2 


(33) 


with 


/*7T 

4, =p+slZ(l)=p+s d(j> cos((j>) R(4 >), 

J — 7T 

A 2 =p+slZ(2)=p+s f dtp cos(20) R(<p). 

J — TV 


The analytical prediction for (x 4 ) n via Eq. (33) is in 
agreement with the simulation results as shown in Fig. 8. 
It can be seen that (x 4 ) n contains similar information as 
the MSD concerning the anomalous diffusive motion of 
the persistent walker. 

The cumulants are often used in the statistical analy¬ 
sis as an alternative to the moments of the distribution. 
In general, the following relations hold between the n- 
th cumulant K n and the moments (shown up to the 4th 
cumulant): 


Hi = (x), 

«2 = ( x 2 )-(x } 2 , 

k 3 = (x 3 )—3(x 2 )(x)+2(x) 3 , 

Ki = (x 4 ‘)—4(x 3 )(x)—3(x 2 ) 2 +12(x 2 )(x} 2 —6(x) 4 . 


(34) 



FIG. 8. (color online). Time evolution of the fourth moment 
of displacement for different values of A, p, and 7 Z. The solid 
lines correspond to analytical predictions and the symbols 
denote simulation results. 


If the walker starts from the origin with the isotropic 
initial condition, the odd moments equal zero and the 
cumulant-moment relations reduce to 


Kl = 0, 

k 2 = {x 2 ), 

k 3 = 0, 

Ki = (x 4 : )—3(x 2 ) 2 . 


(35) 


In the case of an ordinary random walk with 
T(£)=S(£—L), we have p=4,=4 2 =0, (x 2 ) / L 2 =^n, and 
(x 4 )/L 4 =|n 2 —|n, which lead to K 4 /L 4 =— |n (see Fig. 9 
for comparison with simulation). Thus, from Eqs. (21), 
(33), and (35) one can calculate the cumulants, from 
which other useful quantities such as the skewness /3i 
and kurtosis @2 measures can be obtained as 


fa 


fa 


K4_ = ( X 4 ) _ 
K 2 {:X 2 ) 2 


(36) 


In Fig. 10, the time evolution of kurtosis is shown for dif¬ 
ferent turning-angle distributions. For a simple random 
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FIG. 9. The 4th cumulant of x (in units of (£ 4 )) in terms of 
n for a random walk with p=4,=4 2 =0. Inset: The corre¬ 
sponding 4th cumulant of the net displacement r vs. n. 
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1 10 10 2 io 3 


FIG. 10. The kurtosis measure /?2 in terms of the step number 
n for different turning-angle distributions. The lines corre¬ 
spond to analytical predictions via Eq. (36) and the symbols 
denote simulation results. 

walk, /?2 decreases as —3/2 n. Moreover, in anomalous 
diffusive cases, 02 asymptotically converges to zero since 
the long-term behavior is diffusion. 

It is notable that the higher moments are influenced 
by the details of the shape of the turning-angle distri¬ 
bution R{(j>). While the MSD depends only on 1Z=1Z(1) 


(i.e. (cos</)), (x 4 ) is a function of both 1Z(1) and TZ(2) 
(i.e. (cos (j)) and (cos2 cj))). Thus, looking at the behavior 
of the higher moments would reveal the underlying dif¬ 
ferences between turning-angle distributions, which are 
not visible from the MSD results. For example, 1Z equals 
to zero for the three different distributions i?(</>) intro¬ 
duced in Fig. 6(a), thus, their MSD is the same. However, 
their 1Z(2) is 0, —1, and 1 for Ri(q i), -R 2 (</»), and 7? 3 ((ji>), 
respectively. Therefore, one obtains different analytical 
expressions for their higher moments such as (x 4 ). 

Besides the components of the displacement, the net 
distance r of the walker from the origin is also a 
quantity of interest. For a persistent walk with an 
arbitrary turning-angle distribution, so far, there has 
been no exact closed-form expression for (r). For 
approximate expressions (r)^A(r 2 )(l—| ) (short- 

time [13]) and {r)~\\J'x{r 2 ) (asymptotic [14]), one 
deals with the calculation of (r 2 ) and (r 4 ). The 
second moment (r 2 ) can be obtained from Eq. (21) 
since (r 2 )=(x 2 +y 2 )=(x 2 )+(y 2 ). The fourth moment 
reads (r 4 )=((x 2 +y 2 ) 2 )=(x 4 )-|-2(x 2 ?/ 2 )-|-(y 4 ). I 11 general, 
( x 2 y 2 )^(x 2 )(y 2 ) [see Fig. 11(a)]. In order to calculate 
( x 2 y 2 ), one can start from Eq. (3) and follow the analyt¬ 
ical procedure as explained for arbitrary moments (x l ). 
For example, a simple random walk with p=A 1 =A 2 =0 
and A=1 leads to (x 2 y 2 )=|n 2 —vs. (x 2 )(y 2 )=jn 2 . Fi¬ 
nally, the analytical form of (r 4 ) for a walker with con¬ 
stant step size L and isotropic initial conditions is ob¬ 
tained as 


= ~ t L* 


(3 A 


1 +l)n 2 +(A 1 +l)(n—l)n •^iA(^ 2 (^ 1 + j 4 2 ) (A A ) ~KA A) (A A)) 


-4 

+4 

+4 

+8 


i-A (i-A)(i-A)(A-A) 2 

A, (1 -A™) n+A 2 (3 n-n 2 ) -2 A 2 (H/+l) +A 1 (n-1) + 0.5H, (H 1 +l)(n 2 -n-2) 

-2A n + 2 n+A n i +2A\-A\ (n 2 +3n-9) -H 2 (n-2)+^ (n-2)n 2(A 1 +l)A 2 (n-l) 

(TAdJ 5 + (1 - Al )(l-A 2 ) 

A ( A ( 2 A (Hi +2)+2n-l) -2n+5) - (a (2H^+4A 1 +3) +l) A™ 2(A+1) (H 2 -H™ +1 ) 

(IA4J 4 (1-H 1 )(1-H 2 > 

( 2 A - A(A+i)) (A-A +1 ) n -A( n_1 )( 2 A-A(A+ 1 )) 

(i-A) 3 (A-A) 

4AA(A+i)((i~A)( n— 2 )— (A~A _1 )) 2 ( 2 A-A(A+i))a +1 +AA (A-A) 

(i-A) 3 (A-i) +4 (i-A) 2 (A-A) 

8 A( 2 A-A(A+!)) (A-H/+ 2 ) 4AA(A+A) (A-A -1 )" 


(i-A) 4 (A-A) 


(i — A) 2 (A — A) 2 


(37) 


formalism allows for obtaining (r 4 ) in the more general 
case of p^O and even AyM (i.e. variable step lengths) and 
anisotropic initial conditions. One can similarly calculate 


which is confirmed by the simulation data, as shown in 
Fig. 11(b). When p is set to zero, Eq. (37) reduces to 
the expression recently proposed in [18], even though our 
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FIG. 11. (color online), (a) Comparison between ( x 2 y 2 ) (solid lines) and ( x 2 ){y 2 ) (dashed lines) (both presented in units of 
(£ 4 )), obtained from the simulations where the short-time motion is sub (72.= — 0.9) , normal (72=0), or superdiffusion (72=0.9). 
(b) (r 4 ) vs. n from simulations (symbols) or via Eq. (37) (solid lines). 


the cumulants and relative cumulants such as the kurtosis 
for the net displacement r. For example, one finds that 
—K 4 grows as n 2 +n in the simple case of p=A 1 =A 2 =0, 
as shown in Fig. 9 (inset). 


VI. PROBABILITY DISTRIBUTIONS OF THE 
NET DISTANCE AND TURNING ANGLE 
AFTER n STEPS 


In this section, we show how the probability densities 
of the position of the random walker and its orientation 
evolve with time. First, we study the probability distri¬ 
bution P(r) of the distance r of the persistent random 
walker from the origin in simulations. For a simple ran¬ 
dom walk in 2D, the shape of the distribution at step n 
approaches 


P(r) 


2 r 
aDn 


(38) 


in the large n limit (D is the diffusion constant and 
a=^{€)/v). However, the anomalous motion of the per¬ 
sistent walker at short times alters the shape of P(r) 
as well as its propagation speed. In figure 12(a), P(r) is 
plotted at different values of n for p =0 and three turning- 
angle distributions with 72 = —0.9,0, and 0.9. From 
Eq. (38) one expects that all normal-diffusion data col¬ 
lapse onto a universal curve when P{r)-y/n is plotted ver¬ 
sus r/y/n (see the inset). However, the distributions of 
sub and superdiffusion do not follow such a master curve 
at short times. Indeed, P{r) is narrower and the peak 
shifts to the left (right) for subdiffusion (superdiffusion) 
[see Fig. 12(b) (left)]. In the extreme limit of localization 
or ballistic motion, P(r) will be a delta function at r=0 or 
r=n (£), respectively. A similar comparison at long times 
reveals that P(r) broadens slower (faster) than a simple 
random walk in the case of subdiffusion (superdiffusion). 
The shapes are, however, expected to follow Eq. (38), as 


the asymptotic motion is diffusive with different diffu¬ 
sion coefficients obtained from Eq. (19). When scaled by 
y/Un, one finds 

P{r)-\ / Thi~ —-=e (39) 

ay JJn 

thus, we achieve a data collapse for P(r)-y/Dn vs. r/y/Dn 
in the asymptotic regime of the persistent walks when the 
motility is purely diffusive, as shown in Fig. 12(b) (right). 
For random walks in 3D, Eq. (38) should be replaced with 
P ( r )—i_P 

r v yTr (aDn/2) 3 / 2 e 

Finally, we investigate how the introduced angular cor¬ 
relations between successive steps weaken over time. This 
is reflected e.g. in the evolution of the shape of the prob¬ 
ability f n (a) that the direction of motion after n steps 
makes an angle a with the current direction of motion. 
An ordinary random walk is memoryless meaning that 
the walker gets randomized immediately, thus, f n (at) re¬ 
mains isotropic over the whole range of time. However, 
/ n (a) is expected to exhibit anisotropic shapes at short 
times for persistent walks, with a gradual transition to¬ 
wards isotropic distributions in the limit of large n. Let 
us consider the case p =0 for simplicity, i.e. the walker 
turns to a new direction at each step. The probability 
/ n (a) reads 


L( a )=J --Jd<l>vd<j)n - 1 -cj> n -a), 

(40) 

where <j>i denotes the turning angle at step i. Us¬ 
ing the discrete Fourier transform of the delta function 
8{x—a)=A^ XX-oo one obtains the following 

expression for the walks which are symmetric with re¬ 
spect to the arrival direction 


/„(«) 


1 

27r 


1 + 2 ^ 72”(/c) cos (ka) 
k=l 


(41) 
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n=10 n = 4400 




FIG. 12. (color online), (a) Probability distribution of the distance r (in units of {£)) from the origin at n=60,190,1400, and 
4400 (solid, dashed, dash-dotted, and dotted lines, respectively), separately shown for persistent walks with short-time sub 
(left), normal (middle), and superdiffusive motion (right). Insets: Collapse of P(r) ■ y/n vs. r / y/n. (b) Comparison between 
persistent walks with short-time sub (dashed lines), normal (dash-dotted lines), and superdiffusive motion (dotted lines) at the 
early stages of the walk (left) and after a long time (right). Insets: P(r) ■ yjDn in terms of r/y/Dn, where D is the asymptotic 
diffusion coefficient. The solid lines are obtained from Eq. (39). 




FIG. 13. (a) Evolution of the angular distribution f n (a) of the direction of motion a in the lab frame. The walker initially 
arrives along the +x direction. A comparison is made between motions with short-time sub (72.=—0.9), normal (72=0), and 
superdiffusion (7?.=0.9). The lines are obtained from Eq. (41) and symbols denote simulation results. The gray (dotted) lines 
are guides to eye. (b) Probability / n (±-1) of turning to a perpendicular direction after n steps, from simulations (symbols) or 
via Eq. (41) (dashed lines). 


In the case of a=±7r/2, / n (±7r/2) reflects the chance 
of turning to a perpendicular direction after n steps, 
which can be considered as a measure of the coupling be¬ 
tween longitudinal and perpendicular transport. Figure 
13 shows that Eq. (41) is in agreement with simulation 
results. 


VII. SUMMARY AND OUTLOOK 


A persistent random walk model was introduced to 
study the stochastic motion of self-propelled particles. 
By developing a general master equation formalism and 
a Fourier-Z-transform technique it was shown that ana- 
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lytical exact expressions can be obtained for the time evo¬ 
lution of arbitrary moments of displacement. The combi¬ 
nation of self-propulsion and characteristics of step-size 
and turning-angle distributions lead to a rich transport 
phase diagram at short times. The long-time behavior is, 
however, diffusive since the successive step angles in the 
proposed master equation are only indirectly correlated 
over a few number of steps. This defines a time-scale 
between two arbitrarily chosen steps beyond which the 
steps are practically independent of each other. It will 
be interesting to enhance the correlation range e.g. by in¬ 
troducing (anti-)cross correlations between processivity, 
step sizes, and turning angles. For particular functional 


forms of (anti-)cross correlations, one could even obtain 
a stationary increment for the mean square displacement 
(either sub or superdiffusion) over finite time scales as 
observed for the motion in viscoelastic environments. 
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Appendix A: Algebraic coupled equations for the Taylor expansion coefficients Qi(z, a\m) 


Qo(z,a\m) 


zQo,n=o(a\m ) 

z — (p+s TZ(m)) ’ 


(Al) 


„ „ i„^_ zQi,„ = o(a|m) i 1 nia Qo(z,a\m-l)(p+s1l(m-l)) t 1 ^_ ict Qo(z, a\m+l)(p+sll{m+l)) 

vi \Z^OL\Tn) —^ ( / , -n , ' 2^ 


z—(p+slZ(m)) 2 


z — (p+s TZ(m)) 


z — (p+slZ(m)) 


(A2) 


Q 2 {z,a\m) = 


zQ 2 , n =o(a\m) l Q 0 (z,a\m)(p+sll(m)) 


z—(p+slZ(m)) 2 z — (p+slZ(m)) 

{t} 2 ia Qi(z,a\m-l)(p+sTZ(m~l)) (£) 2 _ iQ Qi( 2 ,a|TO+l)(p+s^(m+l)) 


2 > 


z — (p+s1l(rnj) 


2 ) 


z — (p+s7£(m)) 


1 c 2 ia Qo( z > a l TO ~ 2 ) ( P+s K(m- 2)) 1 ^_ 2ia Qo{z , q|w+2) (p+s ^(m+2)) 

4 2 — (jp+slZ(m)') 4 z — (p+slZ[rn)) 


(A3) 


Q 3 (z,a\m) = 


zQ 3 ,n=o(a\m) 3 {£){£ 2 ) Qi(z,a\m)(p+sll(m)) 


z—(p+slZ(m)) 2 (£ 3 4 ) z — (jp+s1Z(rn)) 

^ 3e_ ia Qo(z,a\m+l)(p+sll(m+l)) + 3 ^ ia Qo(z,a\m-l)(p+sTZ{m-l)) 


z — (p+slZ(m)) 


z — (p+slZ(m)) 


1 3ia Qo{z , a|m+3) (p+s TZ{m+ 3)) l^ 3ia Q 0 (z,a\m-3)(p+sTZ(m-3)) 


z — (p+slZ(m)) 


z — (p+slZ(m)) 


3 _ 2ia (e)(£ 2 ) Qi(z : g\m+2)(p+sll(m+2)) 3 2ia (l)(l 2 ) Qi{z,a\m-2)(p+sll(m-2)) 

4 (£ 3 ) z—(p+slZ(m)) 4 (£ 3 ) z — (p+s'Rirn)') 

3 _ ia (£)(£ 2 ) Q 2 {z,a\m+l)(p+sK{rn+l)) 3 ia (l)(l 2 ) Q 2 (z,a\m-l)(p+sn(m-l)) 

2 (£ 3 ) z—(p+slZ(m)) 2 (£ 3 ) 


z — (p+s1Z(m)) 


(A4) 
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zQ 4 ,n=o{a\m) 3 Qo(z, a\m)(p+sTZ(m)) (£ 1 2 )(£ 2 ) Q 2 (z,a\m)(p+sll(m)) 

4 ’ z—(p+slZ(m)) 8 z — (p+slZ(m)) {£ 4 5 6 ) z — (p+slZ(m)) 

1 e - 2 ia Q l w + 2 ) ( p+s K{m+ 2)) 1 c 2 ia a l m ~ 2 ) (p+ s 7l{m- 2)) 

4 2 — (p+s TZ(m)) 4 2 — (p+s7£(m)) 

+ 1 r ^4i a Qo( z ^ a \ rn +' i ) {P+ s ' R ( m + 4 )) 1 c 4i a Qo(^^l w - 4 )(p+ s ^( m - 4 )) 

16 z — (p+slZ(m)) 16 2 — (p+s lZ(m)) 

3 _ ia (l)(l 3 ) Qi(z,a\rn+l)(p+sK(m+l)) 3 ia (l)(l 3 ) Qi(z,a\m-l)(p+sll(m-l)) 

2 (£ 4 ) z—(p+slZ(m)) 2 (f 4 ) 2 — (p+s H{m)) 

1 - 3 ia (l)(l 3 ) Qi(^,a|m+3)(p+s^(m+3)) 1 3ia (l)(l 3 ) Qi(z,a|m-3)(jH-s'fc(m-3)) 

2 (f? 4 ) 2 — (p+s lZ(m)) 2 (^ 4 ) 2 — (p+s7£(m)) 

3 _ 2ia {£ 2 ){£ 2 ) Q 2 (^,Q|m+ 2 )(p+s^(m+ 2 )) 3 2ja (l 2 )(l 2 ) QzQ, a|m-2)(p+s^(m-2)) 
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Appendix B: The first two moments of displacement in the 2-space 
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